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' The formation of molecular hydrogen in the interstellar medium takes 

I place on the surfaces of dust grains. Hydrogen molecules play a role in gas- 

phase reactions that produce other molecules, some of which serve as coolants 
during gravitational collapse and star formation. Thus, the evaluation of the 



o 

i—^ ' production rate of hydrogen molecules and its dependence on the physical con- 

^ I ditions in the cloud are of great importance. Interstellar dust grains exhibit a 

broad size distribution in which the small grains capture most of the surface 

■ area. Recent studies have shown that the production efficiency strongly de- 

5h ' 

' pends on the grain composition and temperature as well as on its size. In this 

paper we present a formula which provides the total production rate of II2 per 
unit volume in the cloud, taking into account the grain composition and tem- 
perature as well as the grain size distribution. The formula agrees very well 
with the master equation results. It shows that for a physically relevant range 
of grain temperatures, the production rate of II2 is significantly enhanced due 
to their broad size distribution. 
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1 INTRODUCTION 



The chemistry of interstellar clo uds consists of reactions taking place in the gas phase as well 



as on the surfaces of dust grains (|Hartquist and William; 



1995^ . In particular, the formation 



of H2 takes place on grain surfaces while all competing gas phase processes are orders of 



2 A. Lipshtat and, 0. Biham 



magnitude less efficient (|Gould fc Salpeteiill963h . Hydrogen molecules are necessary for the 
initiation of gas-phase reaction networks that give rise to the chemical complexity observed 
in interstellar clouds. Molecules serve as coolants during gravitational collapse, enabling the 
formation of stars. 

In diffuse clouds the H2 formation process takes place on bare grains that consist of 
amorphous silicates and carbon. In dense molecular clouds the grains are coated by ice 
mantles. Upon collision with a grain, an H atom has a probability ^ to become adsorbed 
on the surface. The adsorbed H atom (adatom) resides on the surface for an average time 

(residence time) before it desorbs. In the Langmuir-Hinshelwood mechanism, the ad- 
sorbed H atoms diffuse on the surface of the grain either by thermal hopping or by tun- 



neling. When two H atoms encounter each other on the surface, an I: 



(IWilham: 



1968: 



HoUenbach fc Salpeter 



197(1 



1971 



HoUenbach et al 



9. mo 



e cule may form 



1971 



Smoluchowski 



19811 



1983 



Aronowitz &: Chang 



Sandford fc Mlamandola 



The formation rate -Rtotai (cm ^s ^) of molecular hydrogen 



1985; 



Dulev fc Wilhami 



i98e; 



Pirronello fc Avera 



198 



on dust-grain surfaces is commonly evaluated using the HoUenbach- Werner-Salpeter formula 
1 



R 



total 



:^H(^^H)%rain(0"g)7, 



(1) 



where nn (cm~^) is the number density of hydrogen atoms in the gas, (vu) (cm s~^) is the 
average velocity of hydrogen atoms, rigrain (cm ~^) is the number dens ity of dust grains and 
(cg) (cm^) is the average grain cross section (jHoUenbach et al.l 119711 ). Assuming spherical 
grains, (cr^) = '7r(r^) where r (cm) is the grain radius. The parameter ^ 7 ^ 1 is the 
recombination efficiency, namely, the fraction of atoms hitting the grain which come out in 
molecular form. 

The total grain mass amounts to about 1% of the hydrogen mass in the cloud. The 
number density of dust grains is about lO^^'^ — 10^^^ smaller than the total number density 



of hydrogen atoms. In previous studies, a tvpical grain radius o 



a mass density of pg = 2 (gram cm ^) (jHollenbach et al 



0.17/im was used, assuming 



19711 ). Under these assumptions, 



for gas temperature of 100 K and 7 = 0.3, Eq. ((T)) can be simplified to 



total 



R ■ riun, 



(2) 



where R 



10 (cm^ s ^) is the rate coefficient and n = nn + 2nHo (cm ^) is the total 



density of hydrogen atoms in both atomic and molecular (riHa) form ( HoUenbach et al.lll971 ) 



Eq. Q is commonly used for the evaluation of the hydrogen recombination rate in models 



of interstellar chemistry. In this paper we present a direct calculation of -Rtotah without using 
rate coefficients. 

The grain size distribution ng{r) (cm~^) in interstellar clouds is broad with many small 
grains and few large grains. The surface area is dominated by the small grain s . This distri 



butio n can be approximated by a power law of the form (jMathis et al, 
20011) 



ma 



Weingartner 



ng[r) = or 



(3) 



where ng{r)dr (cm~^) is the density of grains with radii in the range (r, r + dr) in the cloud. 
This distribution is bounded between t he upper cutoff rma^ = 0.25/im and the lower cutoff 
'"min = 5nm. The exponent is a ~ 3.5 (jPraine fc Led 119841 ). In general, the normalization 
constant is given by 
3Mg(4 - a) 



rgK' max ' mm / 



(4) 



where Mg (gram cm ^) is the mass density of grain material in the cloud. The total number 



density, ngrain (cm 



of grains is given by 



n. 



gram 



^\ min max J 



a 



(5) 



Recent experiments have shown that the mobility of H atoms on dust an alogues under in 



terstel la r cond i tions is dorninated by thermal hoppi ng rather than tunneling (jPirronello et al 



Il997al lb. 1999! 



Katz et al. 



1999; 



Roser et al. 



20021 ). As a result, the production efficiency of 



H2 molecules is highly dependent on the grain temperature. The efficiency is very high 
within a narrow window along the temperature axis, with exponential decay on both sides. 
Theoretical studies have shown that for sufficiently small grains, the production efficiency 
i s strongly dependent on the grain s ize, namely it decreases as the grain size is reduced 



( Biham et al, 



2001 



Green et al. 



20011 ). The Hollenbach-Werner-Salpeter formula, in which 



all the surface processes are captured in the parameter 7 cannot account for the temperature 
and grain-size dependence of the production efficiency. 

In this paper we present a formula for the production rate of H2 in interstellar clouds, 
which takes into account explicitly the dependence of the recombination efficiency on the 
grain temperature and size, integrating the size dependence over the entire distribution of 
grain sizes. This formula is in very good agreement with the results obtained from the master 
equation. It applies for a broad range of steady-state and quasi steady-state conditions and 
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enables the evaluation of the production rate of H2 in a closed form without the need for 
any simulation or numerical integration. 

The paper is organized as follows. In Sec. 2 we present the different approaches to the 
calculation of reaction rates on surfaces. The dependence of these rates on the grain size 
and temperature is considered in Sec. 3. In Sec. 4 we introduce a formula that accounts for 
the formation rate of molecular hydrogen per unit volume of the clouds, integrated over the 
entire distribution of grain sizes. The results are discussed in Sec. 5 and summarized in Sec. 
6. 



2 RECOMBINATION RATES ON SINGLE GRAINS 



Calculations of reaction rates of Ho formation and other chemical 



are typicallv done using rate equation models ( 



1^ 



Brown fc Charnle^llQQnl: 



Ha,segawa, et a1 



ickles fc Williams 



m 



1 



Caselli et al 



reactions on dust grains 



1977 



d'Hendecourt et al 



WiHacv Williams 



19931 ). These models consist of coupled ordinary differential equations that provide the time 
derivatives of the populations of the species involved. The reaction rates are obtained by 
numerical integration of these equations. The production of molecular hydrogen on a single 
grain of radius r is described by 

^ = F-WN- 2AN\ (6) 
dt ^ ^ 

where N is the number of H atoms on the grain. The first term on the right hand side is 
the flux F (sec~^) of hydrogen atoms onto the grain surface. The second term describes the 
thermal desorption of H atoms from the surface where 

W = uexpi-Ei/kBT) (7) 

is the desorption rate, u is the attempt rate (standardly taken as lO^^sec"^), Ei is the 
activation energy barrier for desorption and T (K) is the grain temperature. The third term 
in Eq. (jUj) accounts for the depletion in the number of adsorbed atoms due to formation of 
molecules. The parameter A = a/S is the rate at which H atoms scan the entire surface of 
the grain, where 

a = u exp^—Eo/ k^T) (8) 

is the hopping rate between adjacent adsorption sites on the grain, Eq is the energy barrier 
for hopping, S = Anr'^s is the number of adsorption sites on the grain surface and s (cm~^) 
is their density. The formation rate, -Rgrain of H2 on a single grain is given by -Rgrain = AN"^. 



Under steady state condit ions one can obtain an exact solution for -Rgrain in terms of F, W 



and A ( Biham et al, 



mi). 



Katz et al. 



( 19991 ) for the amorphous 



Following the experimental results reported by 
carbon sample, we use t he parameters: 5 = 5- 1 0^^ (cm~^), Eq = 44.0 (meV) and Ei = 
56.7 (meV). It is found jRiham TJpshtaJknnA that the process of molecular hydrogen 
formation is efficient only within a narrow window of grain temperatures, Tq < T < Ti, 
where 

Eq 



To-- 
and 
Ti = 



kB HuS/F) 
1E\ — Eq 



(9) 



(10) 



fcB ln(z/5/F) ■ 

At temperatures below Tq the mobility of H atoms on the surface is very low, sharply 
reducing the production rate. At temperatures above Ti most atoms quickly desorb before 
they encounter each other and form molecules. 

Rate equations are an ideal tool for the simulation of surface reactions, due to their 
simplicity and high computational efficiency. In particular, they account correctly for the 
temperature dependence of the reaction rates. However, in the limit of small grains under 
low flux they become unsuitable. This is because they ig nore the fluctuations as well as the 



discrete nature of the populations of atoms on the grain (jCharnlev et al 



1998; 



Shalabiea et al. 



1998; 



Stantcheva et al. 



1997 



Caselli et al 



20011 1. For example, as the number of H atoms 



on a grain fluctuates in the range of 0, 1 or 2, the H2 formation rate cannot be obtained from 
the average number alone. This can be easily understood, since the recombination process 
requires at least two H atoms simultaneously on the surface. 

Recently, a master equation approach was proposed, that takes into account both the 
discrete nature of the population of H atoms as well as the fl uctuations, and is thus suitable 



for th e simulation of H2 formation on interstellar dust grains (iBiham et al.ll2001 



Green et a1 



20011 ). Its dynamical variables are the probabilities P{N) that there are N atoms on the 



grain. The time derivatives P{N), N = 0,1,2, .. ., are expressed in terms of the adsorption, 
reaction and desorption terms according to 



P{N) = F[P{N -1) - P{N)] 

+ [(iV + 1)P(A^ + 1) - A^P(A^)] 
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+ A[{N + 2){N+l)P{N + 2)-N{N-l)P{N)] 



(11) 



where = 0, 1,2, . . .. Note that the equations for P(0) and -P(l) do not include all the 
terms, because at least one H atom is required for desorption to occur and at least two for 
recombination. Direct integration of the master equation provides the time evolution of the 
probabilities P{N). The recombination rate per grain 

Rgrain = A [{N') - {N)) (12) 

can then be calculated, where (A^'^), k = 1,2,..., is the fc'th moment of the distribution 



P(N). Under steady sta te conditions the recombination rate is given by (jGreen et al 



Biham fc Lipshtat 



2m 



"^^^"2V._,(2;^)' ^''^ 

where Ix{y) is the modified Bessel function. 

Using suitable summations over the master equation, one obtains the moment equations, 
namely a set of coupled differential equa tions for the time derivat ives of the moments (A^*^), 
= 1,2, of the distribution P{N) Ihipshtat fc Biham! 12OO3I ). The recombination rate 
on grains of a given size depends only on the first two moments, (A^) and (A^^). The time 
derivatives of these moments are given by 



^^^^ - F+{-W + 2A){N) -2A{N^) 



dt 



d(N'^) 

^ ' - F+ (2F + Vr-4A)(A^) + (8A-2Vr)(A^2) -4A(A^^). (14) 



dt 



Using suitable cutoffs in the master equation one can replace the third moment (A^^) by 
(A^^) = 3{N'^) — 2{N), and obtain a closed set of two coupled differential equations for the 
first two moments. Under steady state conditions Eq. ()14j) takes the form of two coupled 
linear algebraic equations for (A^) and (A^^). The solution of these equations is given by 
(A^) = F{A + W)/{2AF + WA + W^) and (A^^^ = F{F + A + W)/{2AF + WA + W^). 
Inserting these moments into Eq. we obtain 



3 THE LIMITS OF SMALL AND LARGE GRAINS 

Expressing the production rate Rgrain of Eq. (fTHjl in terms of the size-independent hopping 
rate a = AS (sec~^) and the incoming flux / = F/S (monolayer sec~^) we obtain 



if) + + # 



The ratio a/W is the number (up to a lo garithmic correction) of sites that the atom visits 



before it desorbs (MontroU fc Weiss 



19651 ). The ratio W/f is roughly the average number of 
vacant sites in the vicinity of an H atom. Using these parameters we can identify two regimes 
of grain sizes. For grains that are sufficiently large such that S > mm{a/W, W/f), the size 
dependence of the production rate is given by -Rgrain ~ ~ and the production efficiency 
is independent of the grain size. This is the range in which rate equations apply. In the 
limit of small grains, namely when S < min{a/W, W/f), the production efficiency depends 
linearly on S, and -Rgrain ~ 5*^ ~ r^. This is the regime in which the grain size becomes 



the s mallest relevant length-scale in the system and rate equations fail (jBiham fc Lipshtat 



mm. 

To illustrate the implications of these results, consider a given mass M of dust matter. 
Dividing the mass M to spherical grains of radius r, the number of grains is given by 
Rgrain ~ . In an interstellar cloud, the r dependence of the production rate of H2 due to 
such ensemble differs in the two limits discussed above. For large grains the total production 
rate is Rm = ng^^^i^Rgrain ~ f^^^- Reducing the grain size increases their total surface area, 
and as a result the production rate goes up. In the limit of small grains the efficiency 
decreases as the grain size is reduced and thus Rm ~ r. This indicates that there is an 
intermediate grain size for which the production rate is maximal. In Fig. Qwe present the 
total production rate Rm (sec~^) for a total dust mass of M = 1 (gram) which is divided into 
spherical grains of dradius r. Clearly, there is an optimal grain size for which the production 
rate is maximal. This optimal radius depends on the grain composition and temperature 
and is given by 



opt 



'So^t/^ns (17) 



where Sopt = iW/a + 2f /W) ^. For typical interstellar conditions this optimal grain radius 
is of the order of tens of nm, with several thousand adsorption sites. 
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4 THE OVERALL PRODUCTION RATE OF MOLECULAR HYDROGEN 

As discussed above, interstellar grains exhibit a broad size distribution of the form ng{r) ~ 
r~". For a > 2 the surface area is dominated by the small grains, while for a > 3 even the 
grain mass is dominated by the small grains. Therefore, these small grains are expected to 
contribute significantly to H2 formation. Taking this size distribution into account we find 
that the production rate on grains with radii in the range (r, r + dr) is given by 

Cir^~°'dr : small grains 
large grains 



ng{r)Rg,s.m{r)dr 



"'dr 



(18) 



where 



Cl = 

and 

C2 = 



W 



Ansapc 



(19) 



(20) 



+ 2a f 

are obtained directly from Eq. (fTHjl in the limits of S* ^ Sopt and S ^ Sopt, respectively. 
Eq. (fTBj) is expected to be highly accurate in the limit of small grains. For large grains its 
accuracy may be limited because its derivation requires to impose a cutoff on the master 



equation. For large grains the rate equations apply, and therefore one can use 

\ 2 

af 



C2 



1 + 



(21) 



16a \ ' ' V ' "W^ 
which is ob tained directly frorn the exact solution of the rate equation (jH)) under steady state 



conditions teiham et al 



To evaluate the total production rate of molecular hydrogen per unit volume in an 
interstellar cloud, one should take into account the size distribution of dust grains, and 
the dependence of the production efficiency on the grain size and temperature. The total 
production rate is given by 



R 



total 



n„(r)Ra 



{r)dr 



(22) 



where ng{r)dr is the density of grains with radii in the range (r, r + dr), and Rgrain{f) can 
be taken either from the exact result of Eq. (fTSj) or from the approximate result of Eq. (fTBjl . 
Assuming a power law distribution of grain sizes, given by Eq. (jH)), and the result of of Eq. 
(I16|l one can divide the range of integration into two parts, namely the domain of small 
grains, rmin < r < Topt, and the domain of large grains, Topt < r < Tmax- Integrating each 
one of them separately, for a 7^ 3, we obtain 



9 

-Rtotal = g _ ^ ('"opt ~ ^min ) ^ g _ ^ (^max ~ ^opt ) 5 (23) 

In the special case of a = 3, the total production rate is 

^total = Y (ropt - '^min) + ^2 In (r^ax/ ^opt) • (24) 

In Fig. 121 we present the total production rate obtained from Eq. (PH)) for flux / = 
6.88 X 10"^ (ML sec~^), which corresponds to gas temperature of 90 K and hydrogen density 
nn = 10 (cm~'^). The results of this formula are shown for the case in which the parameter 
C2 is taken from Eq. (j20p (o) and for the case in which it is taken from Eq. (j2H) (+). Both 
results are in very good agreement with those obtained from the master equation (solid 
line). The rate equation results (dashed line) are found to significantly over-estimate the 
total production for most of the temperature range. The results obtained (using the master 
equation) under the assumption that the same total grain mass is divided into grains of 
identical size (0.17 /xm) are also shown (dotted line). For the range of grain temperatures in 
which the H2 production is efficient, this assumption leads to an under-estimate of the total 
production rate by about an order of magnitude. This indicates that within this temperature 
range the broad size distribution of the grains tends to enhance the production rate due to 
the increase in grain-surface area. 



5 DISCUSSION 

The formula introduced in this paper [Eq. (jSHI)] enables us to calculate the formation rate of 
molecular hydrogen in interstellar clouds for a broad range of conditions without the need to 
perform computer simulations or numerical integration. It applies both in the case of diffuse 
clouds in which the H2 formation takes place on bare silicate and carbon grains and for 
molecular clouds in which the grains are coated by ice mantles. The formula applies directly 
for grain size distributions that can be approximated by a power-law of the form of Eq. Q. 
For other distribution functions for which the integration cannot be done analytically one 
can solve Eq. ()22|) by numerical integration. 

The relaxation time for II2 formation processes on dust grains can be obtained from Eq. 
(jni). It is given by r = (ly^ + 8ylF)~^/^. Under typical conditions, this time is extremely 
short compared to the time scales of interstellar clouds. Therefore, the formula can be used 
even in the context of evolving clouds. The time dependence can be taken into account by 
updating the input parameters to the formula as the physical conditions in the cloud vary. 
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At very low grain temperatures below Tq, the mobility of H atoms on the surface is 
suppressed and the rate of H2 formation is sharply reduced. In this limit, H atoms keep 
accumulating on the surface, the relaxation time r increases and the behavior may no longer 
be well described by steady state calculations. In this limit the coverage is relatively high 
and thus the rate equations are valid. However, the Langmuir rejection of atoms that hit 
the grain in the vicinity of already adsorbed atoms should be included. This is done by 
replacing the flux term F by F{1 — N/ S). The formula introduced in this paper, Eq. ((221); 
valid only for grain temperatures higher than Tq. To extend its validity below Tq one should 
incorporate the Langmuir rejection term into it, by replacing all the appearances of W in 
Eqs. (Uni) - (EH) by + /. 

In Eq. ()2H|1 it is assumed that all the grains are at the same temperature. This assumption 
does not apply in photon dominated regions (PDR's), where small grains exhibit large 
temperature fluctuations. In general, the formation of molecular hydrogen in PDR's is not 
well understood. The typical grain temperatures in these regions is significantly higher than 
the temperature range in which the formation of H2 on olivine, carbon and ice is found to 
be efficient. Several approaches have been considered in attempt to explain this puzzle. One 
approa ch is based on the assuniption t hat chemisorbed H atoms play a role in the formation 



of H2 (^Cazaux fc Tielens 



20041 ). Other explanations may ivolve porous grains that 



capture H atoms more strongly within pores or with the temperature fluctuations of small 
grains. 



6 SUMMARY 

We have presented a formula for the evaluation of the formation rate of molecular hydrogen 
in interstellar clouds, taking into account the composition and temperature of the grains as 
well as their size distribution. The formula is found to be in excellent agreement with the 
results of the master equation. It enables us to obtain reliable results for the production rates 
without the need to perform computer simulations or numerical integrations. It is found that 
for an astrophysically relevant range of grain temperatures the broad size distribution of the 
grains provides an enhancement of the H2 formation rate compared to previous calculations. 
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Figure 1. The total production rate Rm (sec ^) vs. grain radius r for an ensemble of identical, spherical grains, with total 
grain mass of M = 1 (gram). The density of the grain material is p = 2 (gram cm"^). The grain temperature is T = 16 (solid 
line) and T = 19 (dashed). The flux is 3.44 x 10"^ (ML s-^). The optimal radius is ropt = 5.6 x lO"'^ (cm) for T = 16K and 
T-opt = 1.9 X 10-6 (cm) for T = 19K. 
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Figure 2. The total production rate -Rtotal (cm '^sec ^) of H2 vs. grain temperature as calculated by the formula (Eg. 1231 
introduced in this paper, where C2 is given by Eq. ^ (o) or by Eq. ^ {+). Both results are in very good agreement with 
the master equation (solid line), although the latter one is slightly better. The rate equation (dashed line) is found to over- 
estimate the production rate. These results are compared to those obtained if the same mass of grains is divided into identical, 
spherical grains of 'typical' radium (r = 0.17 n m) (dotted line). It is found that for an astrophysically relevant range of high 
recombination efficiency, the broad size distribution significantly enhances the production rate. 



